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Abstract. In order to take full advantage of the long time series collected by the GONG and 
MDI helioseismic projects, we present here an adaptation of the rotation-corrected m-averaged 
spectrum technique in order to observe low radial-order solar p modes. Modeled profiles of 
the solar rotation demonstrated the potential advantage of such a technique [DOE]. Here we 
develop a new analysis procedure which finds the best estimates of the shift of each m of a given 
(n, I) multiplet, commonly expressed as an expansion in a set of orthogonal polynomials, which 
yield the narrowest mode in the m-averaged spectrum. We apply the technique to the GONG 
data for modes with 1 < £ < 25 and show that it allows us to measure lower-frequency modes 
than with classic peak-fitting analysis of the individual-m spectra. 



1. Introduction 

In the search for long-lived, low-order p modes, the classic methods of peak-fitting the individual- 
m spectra fail because of the decrease in mode amplitudes and increase in solar convective 
motions (solar noise) as the frequency decreases, reducing the signal-to- noise ratio (hereafter 
SNR) of those modes. Instead, various pattern-recognition techniques have been developed in 
an effort to reveal the presence of modes in the low-frequency range. In the case of spatially- 
resolved helioseismic data (for example observations collected by the space-based SOHO/MDI 
instrument and the ground-based multi-site GONG network), m-averaged spectra appeared to 
be a powerful tool. They were employed early in the development of helioseismology |lj, but 
have been replaced by fitting the m-spectra individually as the quality of the data improved. 
Years later, in order to fully take advantage of the long-duration helioseismic MDI and GONG 
instruments, m-averaged spectra corrected by the modeled solar rotation were used to detect new 
low radial-order p modes and to put some upper limits on the detectability of the g modes [21 E]. 
These authors demonstrated the potential advantage of such rotation-corrected m-averaged 
spectra. We present here an adaptation of the m-averaged spectrum technique in which the 
m-dependent shift parameters are determined. 

2. Observation of long-lived, low signal-to-noise-ratio modes 

The accurate observation of low-frequency modes is a difficult task because of two main reasons. 
First, the amplitude of the acoustic modes decreases as the mode inertia increases as the 
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Figure 1. Figures-of-merit (FOM) of the oscillation multiplet t = 5, n = 6 showing the fitted 
peak FWHMs of the m-averaged spectrum as a function of the scanned values of a\, a^, and a§. 



frequency decreases, while the "solar noise" from incoherent, convective motions increases; thus 
the SNR of those modes is progressively reduced. Second, these low-frequency p modes have 
very long lifetimes, as much as several years. These two effects together render the observation 
of such low-amplitude and long-lived modes a difficult task to achieve. However, thanks to 
several long-duration space- and ground-based instruments, more than 11 years of high-quality 
helioseismic data are available today. The observation of low radial-order solar p modes is 
becoming possible if we can somehow improve the SNR of these modes. 

The classic peak-fitting analysis consists of fitting the 21 + 1 individual-m spectra of a given 
multiplet (n, £), either individually or simultaneously. Such fitting methods fail to obtain reliable 
estimates of the mode parameters when the SNR is very low. Instead, for a given multiplet (n, £), 
the mean of 21 + 1 low-SNR spectra can result in an average spectrum with a SNR 3> 1 once 
the individual-m spectra are corrected for the rotation- and structure-induced shifts. A model 
of the solar rotation can be used [H El E] or, as described in this paper, we can also determine 
the shifts as we are searching for the low-frequency modes. 

3. Determination of the splitting coefficients and figures-of-merit 

The m-averaged spectrum is obtained by finding the best estimates of the splitting coefficients, 
commonly called a-coefficients, which yield the narrowest peak in the averaged spectrum. The 
a-coefficients are individually estimated through an iterative process. Then, for a given mode 
(n,£), the frequency shift bv n i m is parameterized by a set of coefficients, as: 

Svntm = ai{n,i)vf\m), (1) 
i=i 

where ai(n,£) are the splitting coefficients, and v\ corresponds to the Clebsch-Gordan 
polynomial expansion [3]. 

For a particular order i of the coefficients cij , a range of values is scanned around the current 
value, while the other ev^s are kept fixed. For each scanned value of aj, the individual-m 
spectra are shifted by the corresponding Clebsch-Gordan polynomials, and the mean of these 
2£+l shifted spectra is taken. The mean of 2£ + 1 independent power densities, which has a x 2 
with more than 2 d.o.f. statistics, can be correctly fitted with a Maximum-Likelihood Estimator 
(MLE) minimization code developed for spectra following a x 2 with 2 d.o.f. statistics if the 
SNR of each m-spectrum is comparable [5]. The mean power spectrum is then MLE fitted by 
a Lorenztian profile and its Full-Width-at-Half-Maximum (FWHM) determined. When all of 
the scans are performed, the corresponding FWHMs can be seen as a figure-of-merit (FOM). 
As shown on Fig. [H the FOMs of the a^s clearly show well-defined minima, the peak profile 
getting narrower as the a-coefficients converge to their best estimates. The uncertainties (cr a .) 
are defined as the formal uncertainties of the fits to the minima. 
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Figure 2. Example of m-averaged power spectra {left panels) before (top) and after correcting 
for the shift coefficients (bottom) in the case of the mode i = 12, n = 4, observed in 2088 
days of GONG data. The estimated central frequencies and a-coefficients are indicated. The 
corresponding m — v diagrams are also shown (right panels). The crosses indicate the position 
of the corresponding theoretical central frequency calculated from Christensen-Dalsgaard's 
model S. 



This iterative procedure is performed until the difference between two iterations in each of 
the computed coefficients aj falls below a given threshold (1CT 2 nHz). A more realistic and 
detailed description of the different effects responsible for lifting the degeneracy of the mode 
frequency into 21 + 1 multiplets needs the use of more splitting coefficients aj in Eq. [T] than we 
have shown here, and is currently underway. 

Finally, the mode parameters (central frequency, linewidth, . . .) are extracted by fitting the 
m-averaged spectrum once the individual-m spectra are corrected by the best estimates of the 
splitting coefficients. However, the formal uncertainties must be normalized by \j2t + 1 [5]. But 
this a-posteriori error normalization is correct only if the 2^ + 1 spectra of a given (n, t) mode 
have the same SNR. The SNRs of the m-spectra inside a multiplet are not the same in our data, 
and the uncertainties of the mode parameters have to be taken as a lower limit only. 

Figure [2] illustrates the advantage of using the m-averaged spectrum technique in the case of 
the mode I = 12, n = 4 at » 1187 /iHz for 2088 days of GONG datail, where the m-averaged 
spectra before and after the correction of the m-spectra for the splitting coefficients are presented. 
The fitted central frequency (i.e. an) and a-coefficients ai,3,5 (determined as on Fig. [IJ) of the 
mode £ = 12, n = 4 are also given on Fig. [2J The corresponding m — v diagram before correction 
(upper-right panel on Fig. [2]) does not show any high SNR structure, while after correction 
(lower-right panel on Fig. [2]), one can see that the individual-m spectra line up. The whole set 
of estimated shifts a\, a%, and as of the low-frequency modes with 1 < I < 25 and measured in 
the 2088-day GONG set are shown on Fig. El 

4. Comparison with other measurements 

In the routine analysis, the GONG and MDI projects use two separate peak-finding pipelines 
to extract the mode parameters. Those pipelines, developed in the early 1990s and mostly 



1 Same observational timespan as used by [8]. 
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Figure 3. Shift coefficients (a±, a^, and 05) of the low radial-order p modes with 1 < i < 25 
measured in 2088 days of GONG data. The a-coefficients are represented as a function of 
frequency v (upper panels), and of vjL, with L = \J 1(1 + 1) (lower panels). The different colors 
correspond to selected ranges of radial order: green dots, modes with n = 1,2; purple, n = 3, 4; 
blue, n = 5,6; and orange, n > 6. 



unchanged since, provide mode parameters on a routine basis. Time series of 108 days are used 
by the GONG project [6], while the MDI project uses 72-day time series [7j. Recently, a new 
and independent peak-finding method of the individual-m spectra, optimized to take advantage 
of the long, spatially-resolved, helioseismic time series available today for both GONG and MDI 
projects, has been developed [8]. 
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Figure 4. i—v diagram of the 
low-frequency modes for t = 
1 — 25 measured in 2088 days 
of GONG data with the m- 
averaged spectrum technique 
(open circles). The crosses 
correspond to the measure- 
ments [S] with the same 2088 
days of GONG and MDI data. 



This new peak-fitting pipeline has been used to extract the low- and medium-degree (t < 25) 
mode parameters on both GONG and MDI observations using one 2088-day long time series, as 
well as using five overlapping segments of 728 days [8j . In order to compare our results obtained 
with the m-averaged spectrum technique, we applied the procedure described in Sec. [3] to the 
same 2088 days of GONG observations. Only modes which passed the so-called HO hypothesis 
statistical test [3] were kept. Figure [H shows the i — v diagram of the low- frequency modes 
measured with the two different analysis. The observed modes with the present method are 
represented by the open circles (GONG data only), contrasted with the crosses (both GONG 
and MDI data) [8]. As seen in Fig. HJ the m-averaged spectrum technique allows us to measure 



a significantly larger number of low-frequency modes in GONG data down to ~ 900 /iHz, and 
to fill in the gaps in the t — v diagram obtained by fitting the individual-m spectra [8]. This 
m-averaged spectrum technique is clearly effective in the low-frequency range where classic peak- 
fitting methods are limited by the low signal-to-noise ratio of each individual-m component. 

5. Conclusion - Perspectives 

We developed a new method to measure low-frequency solar p modes in spatially-resolved 
helioseismic observations. The m-averaged spectrum technique, in which we determine the 
correction for the splitting coefficients (i.e., a-coefficients) , allows us to measure lower frequencies 
in GONG data than with classic peak-fitting methods of the individual-m spectra, which are 
limited by the low signal-to-noise ratio of these low-radial order modes. We are currently in 
the process of generalizing the method with a larger number of a-coefficients, both odd and 
even orders. We are also going to apply this technique to the MDI data. The lower level of 
background noise in the MDI data should allow us to measure even lower low-frequency modes 
than with GONG observations. 

Acknowledgments 

This work utilizes data obtained by the Global Oscillation Network Group (GONG) program, 
managed by the National Solar Observatory, which is operated by AURA, Inc. under a 
cooperative agreement with the National Science Foundation. The data were acquired by 
instruments operated by the Big Bear Solar Observatory, High Altitude Observatory, Learmonth 
Solar Observatory, Udaipur Solar Observatory, Instituto de Astronsica de Canarias, and Cerro 
Tololo Interamerican Observatory. This work has been supported by the NASA SEC GIP grant 
NAG5-11703. 

References 

[1] Brown, T. M. 1985, Nature, 317, 591 

[2] Schou, J. 2004, ESA SP-559: SOHO U Helio- and A stem seismology: Towards a Golden Future, 559, 134 
Schou, J. 2002, ESA SP-508: From Solar Mm to Max: Half a Solar Cycle with SOHO, 508, 99 
Schou, J. 1998, ESA SP-418: Structure and Dynamics of the Interior of the Sun and Sun-like Stars 
SOHO 6/GONG 98, 418, 341 

[3] Appourchaux, T., et al. 2000, ApJ, 538, 401 

[4] Riztwoller, M. H., & Lavely, E. M. 1991, ApJ, 369, 557 

[5] Appourchaux, T. 2003, A&A, 412, 903 

[6] Anderson, E. R., Duvall, T. L., Jr., & Jefferies, S. M. 1990, ApJ, 364, 699 
[7] Schou, J. 1992, Ph.D. Thesis, Aarhus Univ. 
[8] Korzennik, S. G. 2005, ApJ, 626, 585 



